Spin Susceptibility Representation of the Pairing Interaction for the two-dimensional 

Hubbard Model 
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Using numerical dynamic cluster quantum Monte Carlo results, we study a simple approximation 
for the pairing interaction of a two-dimensional Hubbard model with an on-site Coulomb interaction 
U equal to the bandwidth. We find that with an effective temperature dependent coupling U(T) 
and the numerically calculated spin susceptibility x(K ~~ ^Oj ^ ne d-wave pairing interaction is well 
approximated by ^U 2 x(K - K'). 



I. INTRODUCTION 

Numerical calculations have shown that the dominant 
contribution to the d-wave pairing interaction in the Hub- 
bard model comes from a spin S = 1 channel 1,2 . It is 
therefore interesting to determine how well the simple 
RPA form 



°-U\{K-K>) 



(1) 



can describe this interaction. Here x is the spin sus- 
ceptibility, U is an effective coupling strength and K = 
(K.,u! n ) with u> n = (2n + 1)ttT a fermion Matsubara fre- 
quency. Eq. (1) is the form that one expects in weak 
coupling 3-6 . However, we are interested in the case in 
which the Hubbard on-site Coulomb interaction U is 
comparable to the bandwidth, since this is the region 
which is believed to be appropriate for the cuprates 7 . 
This is also the parameter regime in which the pairing is 
the strongest 2 . Here, for the parameters we use through- 
out this paper including a near neighbor hopping t = 1 
and U = 8, the single-particle propagator and density 
of states show clear evidence that one is dealing with a 
doped Mott system. In particular one sees structures as- 
sociated with upper and lower Hubbard bands as well as 
a conduction band 2,8 . In this case, it is unclear whether 
the simple form given by Eq. (1) can adequately describe 
the pairing interaction. We believe that there are two 
reasons that this is possible. First, the leading contri- 
bution of U to the d-wave pairing channel vanishes due 
to the pairing symmetry and secondly, the dressed spin 
susceptibility x an d not the perturbative RPA form will 
be used in Eq. (1). 

To see how well Eq. (1) can describe the pairing inter- 
action, wc first discuss the fitting procedure which will be 
used to determine a temperature dependent effective cou- 
pling U(T). Following this we compare the d-wave eigen- 
value and the K and uj n dependence of the eigenfunc- 
tion of the Bethe-Salpeter equation using the interaction 
given by Eq. (1) with the results obtained using the "ex- 



act" interaction. Throughout this paper, "exact" will re- 
fer to numerical results obtained using a dynamic cluster 
approximation (DCA) Monte Carlo technique 9-12 . We 
will then conclude by discussing the relationship of U to 
the bare interaction. 



II. FITTING THE EFFECTIVE COUPLING U 

For the traditional low T c superconductors, the phonon 
mediated s-wave pairing interaction is characterized by 13 



a 2 F(uj) 



(<5(e P )>p 



(2) 

Here, D v (f\,u)) is the phonon propagator for polarization 
v and {S(e p )} p represents an average over the Fermi sur- 
face weighted by a p-dependent density of states which 
varies as w^ 1 (p). The coupling strength A is given by 



A = 2 / dw 



a 2 F(uj) 



(3) 



and it clearly depends upon the phonon dynamics. How- 
ever, by substituting Eq. (2) into Eq. (3) and making use 
of the usual dispersion relation, A can be conveniently ex- 
pressed in terms of the zero frequency limit of the phonon 
propagator D v (q, uo — 0) 



A 



E 



(|ff 2 . p ,|^(p-p / ,0)^( £p )^( £p Q) pp , 
(<5(e P ))p 



(4) 



At a finite temperature D„(p — p',0) is replaced by 
l/2[£> v (p - p', 0) + D v (p - p', 2*T)]. 

For the Hubbard model we lack a Migdal theorem and 
the pairing interaction is calculated with a dynamic clus- 
ter quantum Monte Carlo algorithm 9-12 . The dynami- 
cal cluster approximation maps the original lattice model 
onto a periodic cluster of size N c sites embedded in a self- 
consistent host. The essential assumption is that short- 
range quantities, such as the self energy and its functional 



2 



derivatives (the irreducible vertex functions) are well rep- 
resented as diagrams constructed from the coarse-grained 
Green's function. For the problem of interest, this is a 
reasonable assumption for systems where the correlations 
that mediate the pairing are short-ranged. To this end, 
the first Brillouin zone is divided into N c cells, with each 
cell represented by its center wave-vector K surrounded 
by N/N c lattice wavevectors labeled by k. The reduc- 
tion of the TV-site lattice problem to an effective N c site 
cluster problem is achieved by coarse-graining the single- 
particle Green's function, i.e. averaging G(K + k) over 
the k within a cell which converges to a cluster Green's 
function G C (K). Consequently, the compact Feynman 
diagrams constructed from G C (K) collapse onto those of 
an effective cluster problem embedded in a host which 
accounts for the fluctuations arising from the hopping of 
electrons between the cluster and the rest of the system. 
The compact cluster quantities are then used to calculate 
the corresponding lattice quantities. 

The pairing interaction is given by the irreducible part 
of the particle-particle vertex 

T PP (K;K') = T PP (K,-K;K',-K') (5) 

with K = (K,uj n ). One can also use the DCA to calcu- 
late the spin susceptibility x(Q, oj n ) 11,12 . In an analogous 
manner to Eq. (4), we now introduce a d-wave coupling 
strength 

l(g(K)r^ en (K, ttT; K', 7rT)ff(K / ))KK' 

(,9 2 (K))k [ ' 

with the even frequency, even momentum part of the 
irreducible particle-particle vertex, 

I*£ ai (K,7rT;K / ,7rT) = ^ ( r w (K, nT; K', ttT) 

+ T pp (K,nT,~-K',irT) 
+ r pp (K,7rT,K',-7rT) 
+ r pp (K,7rT,-K',-7rT)) (7) 

and g(K) — (cosK x — cosK y ). An effective coupling 
U(T) may then be obtained by requiring that the d-wave 
coupling strength in Eq. (6) is the same at a given tem- 
perature for the approximate interaction one obtains by 
replacing T PP (K; K 1 ) by the interaction in Eq. (1). In 
Fig. 1 we show the results for U(T) for two different 
DCA cluster sizes 19 , N c = 4 and N c = 24, for the case 
in which U = 8 and the site filling (n) ~ 0.85. Here one 
sees that U is smaller than U and decreases at lower tem- 
peratures. We will discuss the physics that underlies this 
effect after we explore how well ^U 2 \(K — K') represents 
T pp (K;K'). 
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FIG. 1: The coupling strength U for U = St and a site filling 
(n) = 0.85 calculated for two cluster sizes jV c = 4 and 24. 



III. RESULTS FOR THE PARTICLE-PARTICLE 
BETHE-SALPETER EQUATION 

The leading low temperature eigenvalue of the particle- 
particle Bethe-Salpeter equation 

~ E r even -K; K' , -K') xf W) = 

\MK) (8) 

corresponds to an eigenfunction with d-wave symmetry. 
Here we have coarse-grained the Green's function legs, 
xf{K') = f Ek' G T (K'+k')G ; (-K'-k'), according 
to the DCA assumption. The curves with solid sym- 
bols in Fig. 2 show the d-wave eigenvalue versus T ob- 
tained from Eq. (8) with the "exact" DCA interaction 
T pp . The curves with open symbols show the d-wave 
eigenvalue obtained from Eq. (8) when T pp is replaced 
by ^U 2 x(K — K 1 ). Here we are using DCA results for 
x(K — K 1 ) as well as the single-particle propagator G(k) 
that appears in Eq. (8). One sees that with U(T) deter- 
mined as discussed in Sec. 2, the temperature dependence 
and the size of the d-wave eigenvalue are well accounted 
for by the simple form of the interaction given by Eq. (1). 
For the 4-site cluster, the eigenvalue is larger than that 
for the 24-site cluster because of the absence of pair field 
fluctuations for the 4-site cluster 14 . 

In the following we show results for the 24-site clus- 
ter which was shown in Ref. 14 to be close to convergence 
with respect to the temperature dependence of the pair 
field susceptibility. The momentum dependence of the 
eigenfunction obtained using the approximate form of the 
interaction has the same dominant (cos K x — cos K y ) be- 
havior as the exact DCA result as shown in Fig. 3. Fur- 
thermore, as shown in Fig. 4, the Matsubara frequency 
dependence of the DCA and the approximate interaction 
are remarkably similar. 
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FIG. 2: The d-wave eigenvalue versus T obtained from the 
RPA form, Eq. (1) (open symbols) and from the "exact" DCA 
interaction (solid symbols) for N c = 4 and 24. 
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FIG. 3: The momentum dependence of the d-wave eigenvector 
</3(K,o>„) for ui n = 7rT and T = 0.33 of the approximate 
interaction, Eq. (1), compared to the "exact" DCA interaction 
rfven- <j>(K.,u n ) has been normalized to its value at K = 

M). 



FIG. 4: The frequency dependence of the d-wave eigenvector 
</3(K,a> n ) for K = (tt,0) and T = 0.33 of the approximate 
interaction, Eq. (1), compared to the "exact" DCA interac- 
tion Tfven for two different temperatures. <t5(K, u> n ) has been 
normalized to its value at u>„ — ttT. 



analysis was applied to the heavy fermion superconduc- 
tor UPT3. 15 ' 16 . However, here one would like to have a 
better understanding of U. 




FIG. 5: Decomposition of the irreducible particle-particle ver- 
tex T pp into a fully irreducible vertex Ai rr plus reducible 
particle-hole contributions in the cross channel. Here T 
denotes the full 2-particle vertex and T ph the irreducible 
particle-hole vertex. 



IV. THE EFFECTIVE U(T) 

By fitting U(T) so that the d-wave strength of the 
approximate interaction is equal to that of the "exact" 
DCA interaction, we have found that the approximate 
form, Eq. (1), does an excellent job in describing the d- 
wave eigenvalue and eigenfunction. From a purely phe- 
nomenological point of view, this is an important result. 
It means that to the extent that the Hubbard model 
gives an appropriate description of the cuprates, inelas- 
tic neutron scattering experiments which give 3mx(<7, u>) 
provide a way of determining the momentum and fre- 
quency structure of the pairing interaction. A similar 



As previously discussed 1,2 , the irreducible particle- 
particle vertex can be decomposed into a fully irreducible 
two-fermion vertex Aj rr plus contributions from particle- 
hole channels, as illustrated in Fig. 5. These particle-hole 
channels can be separated into density (S — 0) and mag- 
netic (S — 1) channels. For the even frequency, even 
momentum part of T pp that is of interest, one has 

^event"^' ) ^irr,even 

(K,K') 

+ ^ d (K,K') + ^ m (K,K') (9) 
with the even frequency, even momentum part of the fully 



4 



irreducible two-particle vertex Ai TT>evea 

$ d/m (K,K') = 
1 



and 



(10) 



T r d f m (K - K', uj n - u n , ■ K', w n , ; K, -u n ) 



+ T r d f m (K + K',w„ - w„/; -K',w B /; -K, -w„) 

-ired 



d/r, 

r^/jK-K'.^+^iK', 



-K, — w. 



+ r^(K + K', Wn + w n , ; -K', -w n ,; -K, 



Here, i£ = (K,oj„) and are the two-particle re- 

ducible contributions to the full 4-point vertex V d i m 
calculated in the S=0 (d) or S=l (m) channel, i.e. 



r„ 



where V 



ph 



d/m {Q-K;K') - T p d n /m (Q;K;K') 

d y m is the irreducible vertex in the correspond- 
ing channel. The center of mass and relative wave vectors 
and frequencies in these channels are labeled by the first, 
second and third arguments respectively. In Fig. 6 the d- 
wave projection, Eq. (6), of T^ en (K,K') along with the 
three contributions on the right hand side of Eq. (9) are 
plotted versus the temperature. As expected, the domi- 
nant contribution comes from the magnetic (5 = 1) term. 
The fully irreducible channel is relatively ineffective, but 
at low temperatures the charge density (S — 0) channel 
acts to reduce the effective d-wave pairing strength. This 
is one of the reasons for the decrease found in U(T). A 
similar effect is found in the simple RPA treatment where 



the effective pairing interaction is V sl 



ever, the RPA result for V c h a 
with the temperature. 



charge • 



How- 



does not vary as strongly 



vertex T ph which depends upon K — K 1 as well as internal 
momenta and Matsubara frequencies. In Fig. 7 we show 
results for 



/ (g(K)$ m (K,?rr;^7rr) g (K'))KK' 
^ (,g(K) Xcvcn (K, nT; K', nT)g(K>)} KK , 



1/2 



(11) 



where Xcven is the even frequency, even momentum part 
of the spin susceptibility x obtained in an analogous man- 
ner to Eq. (7). One can see that this estimate for T ph is 
smaller than U. This is consistent with the results of ear- 
lier Quantum Monte Carlo calculations 17,18 which found 
that the electron-spin fluctuation vertex decreased with 
decreasing temperature for all momentum transfers. 
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FIG. 6: Decomposition of the d-wave pairing interaction Vd 
into d-wave projections of the fully irreducible vertex Ai rr , 
and a magnetic (S = 1) channel 3/20 m and a charge (S — 0) 
channel l/2<j>d- 

In addition, from Fig. 5, one can see that even if one 
were to just consider the magnetic channel, U would be 
a more complicated object. The coupling to the S = 
1 susceptibility is given by an irreducible particle-hole 



FIG. 7: Coupling F ph to the S = 1 susceptibility. 



V. CONCLUSION 

To conclude, this work has shown that the momentum 
and frequency dependence of the d-wave pairing inter- 
action for the Hubbard model in the parameter regime 
which is believed to be appropriate to the cuprates is well 
approximated by the spin susceptibility. The strength of 
the coupling when written in the RPA form ^U 2 (T) re- 
quires a temperature dependent effective U(T) which is 
reduced from the bare U and decreases at lower tempera- 
ture in order to phenomenologically account for the effect 
of the charge channel and vertex corrections. 
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